Ground state lost but degeneracy found: 
The effective thermodynamics of artificial spin ice 
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We analyze the rotational demagnetization of artificial spin ice, a recently realized array of 
nanoscale single-domain ferromagnetic islands. Demagnetization does not anneal this model system 
into its anti-ferromagnetic ground state: the moments have a static disordered configuration similar 
to the frozen state of the spin ice materials. We demonstrate that this athermal system has an 
effective extensive degeneracy and we introduce a formalism that can predict the populations of 
local states in this ice-like system with no adjustable parameters. 

PACS numbers: 75.50.Lk, 75.75.+a 



Disordered states are complex and often do not re- 
veal themselves completely to experiment. In neural net- 
works, structural glasses, economic models and countless 
other systems, disorder is often associated with frustra- 
tion, a competition between interactions, not all of which 
can be satisfied. Even perfectly regular lattices can ex- 
hibit frustration if the interactions on the lattice have 
a fundamental geometrical incompatibility, such as for 
antiferromagnetic interactions around a three-fold loop. 
Such geometrical frustration governs the exotic ground 
states of certain spin systems [l|, Q such as spin ice, 
wherein the spin interactions mimic the frustration of 
proton positions in water ice and so produce a degen- 
erate ground state with an extensive zero-temperature 
entropy 0, i, i, H- Recently, Wang et al. 01 have 
fabricated artificial spin ice: a two-dimensional array of 
elongated single-domain permalloy islands whose shape 
anisotropy defines Ising-like spins arranged along the 
sides of a regular square lattice, as in Figs. [T] and [H 
The island-island interactions in such a lattice can be 
engineered to display frustration and the magnetic state 
of every island can be revealed by scanning probe mi- 
croscopy. Frustration-induced disorder on a regular lat- 
tice, fully resolved experimentally, provides a powerful 
model system. As reported in [7[, demagnetized artifi- 
cial spin ice shows short-range ice-like correlations and 
no long-range correlations, with consistent vertex popu- 
lations across experimental runs. A recent paper [8| ex- 
amined this system in the context of the two-dimensional 
vertex models of spin ice [§, [13, U | ■ 

The six pairwise island-island interactions at a given 
vertex of Fig.[2]cannot be simultaneously satisfied. How- 
ever, due to a difference in the pairwise island interaction 
energies for neighbors at tt and 7r/2 angles to each other 
around a given vertex, this array actually has a single 
ground state, unique up to a global spin flip, that carries 
no macroscopic moment. In |7|], the system was demag- 
netized by first fully polarizing it in a large external field 
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FIG. 1: A representative magnetic force microscope image 
of the experimental system. The single domain character of 
islands is indicated by the division of each island into black 
and whites halves 0]. 



and then gradually decreasing the field in an alternating 
stepwise fashion while rapidly rotating the sample, all at 
room temperature [l^ . Although the resulting magnetic 
states had only a small residual moment and a much 
lower energy than the fully polarized starting configura- 
tion, the unique ground state was never reached, nor even 
closely approached. Here, we analyze experimental data 
on artificial spin ice 0] and demonstrate that the demag- 
netization protocol generated a well-defined (albeit not 
thermally equilibrated) disordered state which restores 
the macroscopic degeneracy, but on a hidden manifold. 
The observed vertex populations can then be accounted 
for as the maximally likely outcome within a stochastic 
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FIG. 2: Top; The sixteen possible vertices of artificial spin ice 
and their multiplicities. Bottom: The demagnetized tiling of 
the ground state (left) and a fully polarized tiling of type-II 
vertices (right). 



model of demagnetization. Since the disordered config- 
urations of island moments have very small net magne- 
tization both locally and globally [7|, long-range mag- 
netic interactions are weak and do not sum coherently. 
The dominant interactions within the array connect the 
nearest-neighbor islands that comprise a given vertex. 
Hence we use a 16- vertex model (ol, [lol. . with the ver- 
tices of Fig. [2l the mutual magnetostatic energy of the 
four islands comprising a given vertex is written as Ej, 
Ell, Eiu or Eiv and is denoted as a "vertex energy". All 
vertex energies are calculated with fully relaxed micro- 
magnetic simulations of isolated vertices. For simplicity, 
we set Ej — 0. For the arrays studied, the vertex energies 
maintain the relative order Ei < £"„ < Em < E„ with a 
typical energy scale £'n—i?i ~ 2 x 10^ K (~ 3 attojoules) 
at a lattice constant of 320 nm [l^. This large energy 
scale suppresses thermal fluctuations at room tempera- 
ture, as does the large energy barrier associated with the 
shape anisotropy of the islands, so the spin configurations 
obtained after demagnetization are static. 

If the fractional populations of the four vertex types are 
denoted rij, n„, nm, riiv, then the specific vertex energy 
E — Eini+Eiinii+Emniii+Eivriiv The population frac- 
tions have been extracted by magnetic force microscopy 
on demagnetized arrays. The measured specific vertex 
energy, plotted in Fig. [31 closely tracks that of a pure 
type-II tiling: E « £"„. Even samples subjected to an 
abbreviated anneal protocol, whose residual magnetiza- 
tion is 60% of the saturation value, obey this energetic 
constraint, as does the original fully polarized sample. 
This striking relationship does not arise from simple ran- 
dom averaging- the same plot also shows the average ver- 
tex energy of a random tiling, wherein each vertex type 
occurs according to its multiplicity. Apparently, the ro- 
tational demagnetization protocol, which begins with an 
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FIG. 3: The specific vertex energy E of demagnetized lattices 
(dots) compared to the calculated energy of a type-II ver- 
tex (-Bii, diamonds) and the specific energy of a purely ran- 
dom tiling (triangles) for a range of array lattice constants. 
Error bars incorporate the measured standard deviations of 
ni, riii, nin, riiv across multiple arrays, plus a 5% uncertainty 
in the micromagnetics. The inset shows the specific vertex 
energy and the random tiling energy, normalized to the pure 
type-II vertex energy at each lattice constant. One attojoule 
is 7.25 X 10* Kelvins 



array polarized into a uniform type-II tiling, does not sig- 
nificantly reduce the vertex energy; rather, it moves the 
configuration about within that energy manifold, reduc- 
ing the total magnetostatic energy only through a reduc- 
tion in the long-range demagnetization fields 
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take this experimerrtal observation as given, and explore 
its consequences Population fractions for different 

demagnetized samples at the same lattice constant are 
very similar, presumably reflecting some well-deflned un- 
derlying state. Individual islands for samples with lattice 
constants of 320 and 560 nm show no significant correla- 
tions in moment orientation between successive demag- 
netization runs, so the moments are not pinned by static 
structural disorder in island shape. All samples analyzed 
below had a residual magnetization Mr < 10%. Only 
at the largest lattice constants, where the island-island 
interactions are weakest, does the manifold break down 
and the specific energy approach the random value, i.e. 
that for equal probabilities for each of the 16 vertices. 

To characterize the demagnetized state of this ather- 
mal system, and to predicts its statistics, we attempt a 
description as the most likely outcome of demagnetiza- 
tion process. During demagnetization, the rotating sam- 
ple is placed in a time varying magnetic field [l^ . At first, 
the field is large enough to polarize all the islands into 
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a type-II configuration. As the field magnitude drops, 
it carves defects inside that background. These defects, 
treated as a non-interacting gas, may appear as types I, 
II, III or IV, but the background within which they hve 
is purely type-II. We introduce p, the density of defected 
vertices; it relates the fractional populations of each ver- 
tex type within the defected population, written as i^i, 
Uii, Viii and i>iv, to the total populations rii, n„, n„i, 



riiu = pvu 

riiv = pVl^ 



(1) 



We adopt a vertex-gas approximation [l6[ , wherein each 
vertex is treated as an independent entity. Thus there 
are 
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ways to tear D — pN defects in the N vertices of a po- 
larized tiling, allocated among the four vertex types with 
distribution Na ~ pv^, a = I, . . . , IV. We model the 
demagnetization protocol as a one-step non-equilibrium 
stochastic process at the vertex level. To obtain the most 
likely outcome of the vertex populations after demagne- 
tization, we maximize the logarithm of the multiplicity 
M of Eq. H 



s = - [pln/7+ (1 - p) ln(l - p)] + per, 

formally an entropy, under an energy constraint E 
The second term of Eq. [3] 
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is formally the entropy of an ensemble of vertices a, each 
with probability Va and given multiplicity. However, if 
taken as the real entropy for the actual vertex population 
TT-Q, or —111 In irii - n„ In in„ - ni„ In in„i — n„ In ^niv, 
and maximized, a would return populations in only mod- 
est agreement with the experimental results for type-II, 
as seen in Fig. [H As already stressed, demagnetization 
is not thermal equilibration. Maximization of the "en- 
tropy" s of Eq. [3] simply returns the most likely outcome 
of the demagnetization. 

Since background vertices all contribute i?„ to the en- 
ergy, the energy constraint is relevant only in the maxi- 
mization of the "inner entropy" tr: £'ii^i+£'„f„+£'nii^iii + 
EivVi-v = Ell. The defect population fractions that result 
are written {t^*}. Next, s is maximized relative to the 
total defect population p with cr evaluated at {i'*}, yield- 
ing 
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FIG. 4: Experimental population fractions of vertex types 
(dots) compared to the values predicted by maximizing the 
entropy on the type-II vertex energy manifold, with (dia- 
monds) and without (triangles) the background contribution. 
Dashed lines give the population fraction of a random tiling. 
Error bars represent statistical counting errors in the popula- 
tions averaged across sample runs. 



The residual magnetization (as a fraction of the satu- 
ration magnetization, which is pure type-II) is bounded 
above by the background fraction: Mr < 1 — p 17]. 

The vertex population fractions obtained by this pro- 
cedure are compared to the experimental results in Fig. 3] 
The agreement is good, with no adjustable parameters. 
This agreement, if taken to verify the model, implies a 
striking feature in the process of rotational demagnetiza- 
tion: in probabilistic terms, the background state stands 
on equal footing with each of the defect states. Con- 
sidering the history of each island's magnet state during 
demagnetization (in the rotating frame), this suggests 
that each vertex has a single opportunity to leave the 
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background state and all histories are equally probable. 

As the lattice constant increases, the energy difference 
between types I and II increases faster than that between 
II and III. Therefore, under the energy constraint, the ra- 
tio n,„/n, increases: the defected sample acquires larger 
multiplicity and p must increase; hence type-Ill increase 
and type-II decreases. The fraction of type-IV is too 
small to make strong statements, although its popula- 
tion increase at larger lattice constants is also predicted 
by theory. 

The predictions are least accurate at large lattice con- 
stant, where island- island interactions are weak and the 
vertex approximation breaks down, since the system ap- 
proaches a limit of ideal independent dipoles. At large 
lattice constants, the system is observed to deviate from 
the manifold of constant vertex energy (see inset of 
Fig. [3]), approaching the random-tiling values expected 
for non-interacting islands. Eqn. [5] is therefore expected 
to not apply in this regime. Nevertheless, it can still 
predict at least the population of type-II vertices, which 
are those most affected by p: assuming that the ran- 
dom tiling fraction holds for the defected vertices (i.e. 
I'll = 1/4), Eqn. [5] imphes p = 16/17, so we obtain 
n„ = l — y|-(-y||~ 0.294 . . . , close to the experimental 
values of 0.305 ± 0.013 for 680 nm and 0.287 ± 0.028 for 
880 nm. These arrays appear to occupy an intermediate 
regime where the energy constraint has broken down, but 
not yet the stochastic hypothesis. 

In summary, we have shown that artificial spin ice, 
despite the absence of thermal fluctuations, can be de- 
scribed within an entropy maximization formalism, re- 
sulting in effective thermodynamic behavior for an ather- 
mal system. A demagnetized state with an extensive 
frozen-in residual entropy (similar to that seen in water 
and spin ice) lives on a sub-manifold of constant vertex 
energy as the most likely outcome of a minimal stochas- 
tic model of the demagnetization process. All nontrivial 
structure in the microstate can be subsumed into the 
residual polarized component, with the remainder dis- 
tributed according a Gibbs-type ensemble. The exper- 
imental results are closely matched with no adjustable 
parameters. The dynamical origin of this behavior re- 
mains to be explained, and may relate to the step size 
chosen experimentally being the maximal decrement in 
absolute field that demagnetizes the array. Behavior at 
smaller step sizes could reveal further details as to the 
full extent of this regime. 
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